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Abstract. 

75 As and 31 P NMR studies are performed in PrCoAsO and NdCoPO respectively. 

The Knight shift data in PrCoAsO indicate the presence of an antiferromagnetic 

interaction between the 4/ moments along the c axis in the ferromagnetic state of Co 3d 

moments. We propose a possible spin structure in this system. The 75 As quadrupolar 

q ■ coupling constant, vq increases continuously with the decrease of temperature and 

is found to vary linearly with the intrinsic spin susceptibility, Ki SO , This indicates 
a possibility of the presence of a coupling between charge density and spin density 
fluctuations. Further, 31 P NMR Knight shift and spin lattice relaxation rate (1/Ti) 
in the paramagnetic state of NdCoPO indicate that the differences between LaCoPO 

\^q • and NdCoPO with SmCoPO are due to the decrement of inter layer separation and 

not due to the moments of 4/ electrons. Nuclear spin lattice relaxation time (Ti) in 

■h- | NdCoPO shows weak anisotropy at 300 K. Using self consistent renormalization (SCR) 

^^ ■ theory of itinerant ferromagnct, it is shown that in the ab plane, the spin fluctuations 

are three dimensional ferromagnetic in nature. From SCR theory the important spin 
fluctuation parameters (To, Ta, -Fi) are evaluated. The similarities and dissimilarities 
of the NMR results in As and P based systems, with different rare earths have also 

S^ ' been discussed. 
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1. Introduction 

The main attraction of the layered FeAs based superconductors to the condensed matter 
physicists is due to their unconventional nature and the presence of several competing 
interactions. These superconductors are grouped in several families pQ. In 1111 and 122 
families, superconductivity can be achieved by Co doping in place of iron [21 13] . Whereas, 
in LaCo2B 2 , which also belongs to 122 family, superconductivity appears upon Y doping 
at the La site jl]. Thus it is emerged that Co is playing different role in different 
systems of the same type of crystal structure. Therefore, the non superconducting 
Co based members in 1111 and 122 family should be explored systematically with 
microscopic tools. This may help to understand the physics behind the emergence 
of superconductivity in non-doped cases. Apart from this fact, the novel and complex 
magnetic properties of the Co based magnetic systems RECoAsO/ RE CoPO (where 
RE is rare-earth ion) belonging to 1111 family, are also of interest because of the 4 /-3d 
interplay [HIE]. Surprisingly, in REFeAsO series no such dominant effect of the interplay 
of 3d-4f interaction was observed in magnetic properties. 

Recently it was shown from magnetization and specific heat data that Sm/NdCoPO 
undergo three magnetic transitions i.e. Tc,co (80 K), the Sm 4 -^-Co M and Nd 4 ^-Co M inter- 
played antiferromagnetic (AFM) transition below 20 K and finally Sm 3+ and Nd 3+ spins 
individual AFM transitions at 5.4 K and below 2 K respectively jSJ. Again in RE CoAsO 
series La, Ce, and Pr show PM to FM transition, Sm, Nd, and Gd show PM — y FM 
—y AFM transition J6J Ej. Furthermore, in SmCoAsO and NdCoAsO a second AFM 
transition only due to rare earth ion as in NdCoPO and SmCoPO was also reported 
[8j [9] . The main difference between P and As based compounds is that the FM transition 
temperature (T c ) increases from La to Ce and for Pr, Nd, Sm it is almost same in 
RE CoAsO series whereas in RE CoPO family T c increases progressively as we go down 
the series from La to Sm. In each case, the lattice volume decreases across the series. 
In general, with the application of chemical or physical pressure, T c decreases due to 
the increment of density of state (DOS) at the Fermi level (magneto- volume effect). 
However, due to the lattice size decrement, the three dimensionality of the magnetic 
interaction may enhance causing an increment of T c and also T^ • These point towards 
the active role of the competing phenomena governing the actual ground state and most 
importantly the role of competing 4/ and 3d moments makes these families a novel and 
complex one in the field of magnetism. 

Our earlier 31 P and 139 La NMR measurements in LaCoPO (quasi 2D fermi surface) 
revealed that the spin fluctuation of 3d electrons in the paramagnetic (PM) state is 
basically two dimensional (2D) in nature with non negligible three dimensional (3D) 
part and it is 3D in the ferromagnetic (FM) state [TOj [11]. In SmCoPO, we have 
observed that the 3<i-spin fluctuations in the ab plane is primarily of 2D FM in nature, 
while along the c-axis, a signature of a weak 2D AFM spin fluctuations superimposed 
on weak FM spin-fluctuations, even in a field of 7 T and far above T/v, was observed. 
The interaction between Sm-4/ and Co-3<i spin fluctuations, has a contribution for the 
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development of weak AFM spin fluctuations along c-axis at a temperature far above the 
AFM transition temperature of the Co-3d spins below their FM transition [12]. The 
moment of Sm 4/ and the separation between the ab planes are less in case of SmCoPO 
than in NdCoPO. So, it would be interesting to study NdCoPO to understand the role of 
the change of the inter-plane separation and the moment of 4/ electrons on its complex 
magnetism. Further the nature and the dimensionality of 3d spin-fluctuations within 
the ab plane in NdCoPO may also be compared with that in SmCoPO and LaCoPO. 

Recent muon-spin rotation and relaxation studies in PrCoAsO indicated a further 
change in the magnetic state below the Co FM transition (Tc~ 75 K) [13]. In contrast, 
in PrCoPO, no such change was reported [H]. Since NMR is a sensitive tool to probe 
the magnetic interactions, we intend to study the magnetic behavior of PrCoAsO and 
NdCoPO using 75 As and 31 P NMR respectively. Moreover, for NdCoPO, till now, no 
results on microscopic measurements are reported. Also the 75 As NMR results are 
reported for the first time in this paper. Finally, the present results are compared with 
those in La based analogue, in order to understand the role of 4/ electrons in As and P 
based series. 

2. Experimental 

For the synthesis of polycrystalline sample of PrCoAsO, at first, ingot of Pr metal (Alfa 
Aesar, 99.9 %) is cut into small pieces and grounded to make tiny pieces. It is then 
mixed with As powder(Alfa Aesar, 99.99 %) and the mixture is kept inside a silica tube 
to prepare PrAs. The whole process is done in an argon atmosphere. The silica tube is 
evacuated and sealed and then carefully fired in a furnace at 550 °C for 5 h and then at 
800 °C for 12 h. The obtained powder of PrAs is mixed with the powders of CoO (Alfa 
Aesar, 99.999 %) in stoichiometric ratio and ground well in argon to avoid oxidation. 
The pelletized mixture of PrAs and CoO is then kept into an evacuated silica tube and 
fired at 1100 °C for 12 h. The powder sample is characterized by using x-ray diffraction 
technique with CuKa radiation at room temperature in a Rigaku X-ray diffractometer. 
The Rietveld analysis of the X-ray powder diffraction data (figure 1) confirmed that 
the sample is crystallized in the tetragonal phase with majority of the peaks indexed 
to the space group P4/nmm with lattice constants a = b = 4.001 A and c = 8.35 A. 
There are some unindexed lines with small intensity which are arising from the impurity 
phases of Praseodymium oxide and CoAs. The procedure of preparation of NdCoPO 
was discussed in ref: 6. 

The 75 As and 31 P NMR measurements are performed on the powder samples of 
PrCoAsO and NdCoPO using a conventional phase- coherent spectrometer (Thamway 
PROT 4103MR) with a 7.0 T superconducting magnet (Bruker). The spectrum is 
recorded by changing the frequency step by step and recording the spin echo intensity 
by applying a 7r/2 — r — n/2 solid echo pulse sequence. The temperature variation study 
is performed in an Oxford continuous flow cryostat equipped with a ITC503 controller. 
The spin lattice relaxation time (7\) is measured using the saturation recovery method, 
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Figure 1. (Color online) Powder X-ray diffraction pattern of PrCoAsO at 300 K with 
Rietveld fit. 



applying a single n/2 pulse. 



3. PrCoAsO 



3.1. 75 As NMR spectra 




Frequency (MHz) 

Figure 2. (Color online) 75 As NMR spectra at 300 K in PrCoAsO. The continuous line 
corresponds to the theoretical fit. The vertical dashed line corresponds to the reference 
position of 75 As nuclei at the field of 7T. Inset (a) shows temperature dependence of 
v Q (MHz) for PrCoAsO, inset (b) shows v Q (MHz) versus V (A 3 ) for LaCoAsO [37], 
CeCoAsO El and PrCoAsO at 300 K. 



Figure 2 shows the 75 As (1=3/2) NMR spectrum at 300 K in PrCoAsO which 
consists of a splitted central line and two satellites. The splitting of the central line is due 
to the combined effect of magnetic anisotropy and second order quadrupolar interaction. 
Assuming the principal axes of the electric field gradient (EFG) and magnetic shift 
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Figure 3. (Color online) 75 As NMR spectra at different temperatures in PrCoAsO. 
The vertical dashed line corresponds to the reference position of 75 As nuclei at the 
field of 7 T. 



tensors are parallel to each other, the spectrum can be simulated as shown by continuous 
line in figure 2. 

Unlike the case of CeCoAsO [5J, 75 As NMR spectrum in PrCoAsO (figure 3) did 
not vanish below the FM transition temperature (Tc ~ 75 K). However, the magnetic 
anisotropy and the line width increase continuously. The increase in line width can be 
understood as due to the effect of the increase of magnetization. 31 P NMR spectrum 
in LaCoPO pU] was also detected below Tc- Whereas, in SmCoPO where Tc=110 K 
in a field of 7 T, the NMR line could not be detected below 130 K, due to the large 
increase in line width and anisotropy of the hyperfine field from below 300 K, originating 
mainly from Sm 4/ spin contribution [12]. This indicates that in PrCoAsO the effect of 
4/ moment on 75 As NMR line shape is weaker than that in SmCoPO. The 75 As NMR 
results in SmCoAsO are not yet reported. From figure 3 it is seen that no noticeable 
change in the 75 As NMR line shape and position in PrCoAsO are discernable below 40 
K, as expected from a change in the magnetic state, observed in /xSR results below 50 
K [13] . A possible reason for this could be the suppression of the change in the local 
magnetic field below 50 K (observed from /iSR study in zero magnetic field) in a field 
of 7 T, used in present NMR experiment. Recently from /iSR, it was reported that the 
results of PrCoPO (T c =48 K) are qualitatively identical to that of LaCoPO, despite 
the presence of magnetic moment of Pr 3+ ion [T4"] . 

In the paramagnetic state of PrCoAsO, 75 As nuclear quadrupolar coupling constant 
(uq) increases continuously below 300 K, (shown in inset (a) of figure 2) whose origin 
will be discussed later. This kind of increment of vq was also observed in 75 As NMR of 
CeCoAsO [5]. In LaCoAsO [15], CeCoAsO [5] and PrCoAsO, the values of v Q at 300 
K are 2.9, 3.6 and 3.9 MHz respectively indicating that the local EFG at the 75 As site 
increases as we go from La to Pr. The distortion of the Co-As tetrahedra is related to 
decrement of the lattice volume from La to Pr as can be seen from the linear variation 
of Uq's (near room temperature) with the unit cell volume (inset (b) of figure 2). 
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3.2. Knight shift and hyperfine field 

By fitting the experimental spectra with the simulated theoretical line we have estimated 
K c , K a i,, Ki SO , K ax and vq. K c , K ab are the Knight shifts corresponding to #=0 
and 8=ir/2, where 8 is the angle between the direction of the external magnetic 
field and z principle axis of the magnetic/electrostatic hyperfine interaction tensor. 
Kiso = 2/3K ab + l/3K c and K ax = 1/3(K C — K ab ). Figure 4 and 5 shows the temperature 
dependence of shift parameters K c , K ab , K iso , K ax in PrCoAsO down to 10 K and for 
comparison, those reported in LaCoAsO [15] and CeCoAsO [5] have been included. 
The measured shift, K = K Q + K(T), where K is the temperature independent 
contribution arising from conduction electron spin susceptibility, orbital susceptibility 
and diamagnetic susceptibility of core electrons. K(T) arises from the temperature 
dependent susceptibility due to Go-3d and Pr-4/ spins, 

K(T) = (H hf /Npi B ) X (T) (1) 

H h f is the total hyperfine field related to the total hyperfine coupling constant A h f by 
Hhf = Ahf/'yhg, N is the Avogadro number and //# is the Bohr magneton. Inset of 
figure 5 shows the linear variation of K iso and K ax with x = M/H. The estimated 
values of H^f and Hffi are 15.35 and -2.5 kOe//i£ respectively. 

Using the structural parameters for PrCoAsO we have calculated the dipolar field 
at the 31 As nuclear site arising from Pr moment with the formula, 

v^ (3rjr fc - r 2 5 jk ) . 
H dip = ^2^ ~ 5 ]J,k = x,y,z (2) 

jjl is the magnetic moment of the Pr ion. The H^L (-0.11 kOe/yUs) is one order of 
magnitude smaller than the experimental value of Hffi. This indicates that the measured 
Hf5 is governed mainly by the anisotropic hyperfine field produced by Pr-4/ and Co-3d 
moments at the 75 As nuclear site. So the observed values of H h f could be related to 
the effect of hybridization of the 4s and 4p-orbitals of As with the 3<i-orbitals of Co and 
4/-orbitals of Pr via the conduction electrons. 

Temperature dependence of K iso in the paramagnetic phase can be well described 
by the Curie- Weiss type behavior, 

C 

if 

as represented by the continuous lines in figure 5 for LaCoAsO, CeCoAsO and PrCoAsO. 
The estimated value of the effective moment, P e // from the Curie- Weiss constant (C) in 
LaCoAsO is 1.64/i B with the Curie- Weiss temperature, #=70 K. The reduced value of 
Co effective moment with respect to the local Co moment (3.87/x.b) indicates that Co 
3d electrons are of itinerant in nature. For CeCoAsO and PrCoAsO these parameters 
are 2/^g (with 8=75 K) and 2.18/ig (with 8=80 K) respectively. The P e ff in case of 
PrCoAsO and CeCoAsO systems is due to the contribution of Co 3d and rare earth Af 
i.e. P e ff = \/(fxfrf) 2 + (fJ-lff) 2 - If we assume that for all other rare earths the P e ff for 
Co 3d is same then it is clear that in the paramagnetic state, there is also a contribution 
to P e ff from 4/ electrons over that of the Co 3d electrons. 



K iso (T) = (fl357/^A*B)7^— 5T ( 3 ) 



PrCoAsO and NdCoPO 





■ 


1 * 


i. 


20- 




* 


- 


16- 






.. 




•»g 




■ 


* 12- 






- 






s„ 


. 


/ 8- 




1. 


- 






B* 


■ 


4- 




*' 


_ 






* •#• » a* a* » 


a* 


40- 




• 








* 


■ 


30- 






- 






* 


■ 


£. 20- 






- 


■£. 




* 


■ 


10- 


t»)9 




. 


0- 


H H 


H *» ft »-» *•*«*» 


j* 





50 100 150 200 250 


300 



Temperature (K) 



Figure 4. K c (%), K a b (%) versus temperature for LaCoAsO (star), CeCoAsO (half 
filled square) and PrCoAsO (filled circle) . Data for LaCoAsO and CeCoAsO have been 
taken from [37j and [5] respectively. 
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Figure 5. Ki SO (%), K ax (%) versus temperature for LaCoAsO (star), CeCoAsO (half 
filled square) and PrCoAsO (filled circle). Data for LaCoAsO and CeCoAsO have been 
taken from [37] and [5] respectively. Inset shows Ki SO , K ax versus M/H for PrCoAsO. 



H™? and Hfij are found to vary linearly with the unit cell volume of LaCoAsO, 
CeCoAsO and PrCoAsO (figure 6). This indicates that the extent of hybridization of 
the electronic orbitals, involved in the hyperfine coupling constant, changes linearly with 
the change in the volume of the unit cell. Further H l h s ? also scales with P e ff (figure 6). 
As the values of the spontaneous magnetization M s at T=0 and the Tc are almost same 
in CeCoAsO and PrCoAsO [7] the enhanced P e // values for CeCoAsO and PrCoAsO 
compared to that in LaCoAsO should be due to the rare earth-4/. Also the P e ff values 
are less than those corresponding to the localized 4/ moment of Ce and Pr. This could 
arise due to the partial delocalization of 4/ electrons, due to a change in the band 
structure, resulting from the decrease in the unit cell volume for substitution of Ce and 
Pr in place of La. 
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Figure 6. (Color online) (a) Wfi? versus V (A 3 ) and (b) H™? (closed symbol) and 
Hff (open symbol) versus P eff (fj, B ) for LaCoAsO [37], CeCoAsO [5] and PrCoAsO. 



The values of K a b for LaCoAsO, CeCoAsO and PrCoAsO are almost same in the 
paramagnetic state but K c is lower in CeCoAsO and PrCoAsO compared to that in 
LaCoAsO which suggests that the local field produced at the 75 As site along c-direction 
due to 4/ electrons, is of opposite sign compared to that of Co-3d. In case of CeCoAsO 
a peak in shift was observed near at 50 K, which was related to the crystal electric field 
(CEF) splitting of Ce 4/ and to the Schottky-type anomaly in specific heat [5]. No such 
peak is observed in PrCoAsO. The present Knight shift data suggest that even at 80 K 
the 75 As nucleus experiences a local magnetic field along c axis due to the 4/ moments, 
because within the ab plane, this field is same for La and Ce based analog, whereas, 
along c axis a drop in magnitude was observed. The inter layer AFM correlation of 
Co 3d or Pr 4/ and intralayer FM correlation of Co 3d moments can be the moment 
configuration. If this happens then in the Co-As layer the local field along c axis can 
be small due to the cancelation of AFM moments. This was also seen in NdCoAsO [9] 
and SmCoPO |12j . Further experiments like neutron elastic scattering have to be done 
to know the exact reason for this. The effect of CEF on the Pr-4/ levels can not be also 
neglected for the decrement of shift compared to the La analogue. 



3.3. Possible origin of the temperature dependence of quadrupolar coupling constant 

In case of PrCoAsO we have seen from 75 As NMR that the Uq increases with the 
decrease of temperature (inset (a) of figure 2). The thermal (phonon) contribution can 
induce an increase in vq at lower temperature. In this case, it follows the relation 
v onon _ u o_ a y3/2 w ]j ere u o j g ^ e uq 8kt T — K. It can easily be seen from the inset 
(a) of figure 2 that the temperature dependence of vq does not follow this power law. 

Most interestingly, we have seen that vq scales with K iso or x (figure 7) . This type 
of scaling was seen in case of MnSi (itinerant ferromagnet) [16]. To explain the linearity, 
a theory based on self consistent renormalization (SCR) approach was developed [T7] . 
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In general, the total electric field gradient (q) in metals can be written as 

/■oo 

q= dr(l + 7(r))g(r) (4) 

Jo 

where q(r) is the quadrupole charge density between r and r + dr and 7(7*) is its 

antishielding factor at r. The electric field gradient q in metals consists of two 

contributions, 

q = qi + q e i(T) (5) 

where qi is the ion core (lattice) electric field gradient and q e \ is the non-cubic part of 
the conduction electron charge distribution. One can easily evaluate the magnitude of 
qi by extrapolating the linear portion of uq versus x or K iso plot to x or K iso = 0. A 
contribution to q e \ in terms of quadrupolar charge susceptibility, which can be influenced 
by spin-fluctuations through the mode-mode coupling between charge and spin density 
fluctuations, was calculated and for which a linear relation between uq and x was found 
|17j . Thus the linear relation between uq and the intrinsic spin susceptibility K iso in 
PrCoAsO indicates a possible signature of the presence of the coupling between charge 
density and spin density fluctuations in PrCoAsO resulting a temperature dependent Uq. 
Since UQ. e i(T) is proportional to x an d not to magnetization, so it should not change 
with magnetic field. It is also necessary to do the NMR experiments in presence of 
different magnetic fields, in order to check that uq is field independent, to discard the 
possibility of the presence of static magneto-elastic coupling which can change q^ or q e i, 
as the magnetization increases with increasing magnetic field. 




X (emu/mole) 



Figure 7. (Color online) vq versus local susceptibility for PrCoAsO, BaFe2As2 [18] 
NaCo 2 4 M and LiFeAs [E]. 



We found that itinerant systems like BaFe 2 As 2 [18], LiFeAs [19] and NaCo 2 4 [20] 
also follow the same scaling relation (figure 7) wherein the temperature dependence of 
uq were also observed. In LiFeAs, the value of 75 As uq for H=0 T (from NQR) is 21.12 
MHz (at 20 K) which is same when NMR has been done at 51.1 MHz [21]. In BaFe 2 AS 2 , 
75 As NMR u Q = 3 MHz, at the resonance frequencies u R = 45 MHz [22] and 48.4MHz 
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j. This field independent value of vq further supports the applicability of the above 
mentioned model to explain the T dependent vq in these systems. 



4. NdCoPO 

4-1. 31 P NMR spectrum, Knight shift and hyperfine field 
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Figure 8. (Color online) 31 P NMR spectrum of NdCoPO at 300 K. Continuous line is 
the total theoretical fit which consists of pure NdCoPO (dashed curve) and impurity 
(dotted curve). The vertical dashed line corresponds to the reference position of 31 P 
nuclei at the field of 7 T. 

Figure 8 shows the typical 31 P NMR spectrum in polycrystalline NdCoPO at 300 
K. The resonance line shape corresponds to a powder pattern of a spin I = 1/2 nucleus 
experiencing an axially symmetric local magnetic field, as expected for NdCoPO having 
tetragonal symmetry. The step in the low-frequency side corresponds to H \\ c (6 = 0°) 
and the maximum in high frequency side corresponds to H Q _L c (6 = 90°). The shift of 
the step with respect to the reference position {vr), corresponds to K c and that of the 
maximum corresponds to K a jj. The continuous line superimposed on the experimental 
line is the calculated spectrum using Gaussian line shape along with the symmetric 
impurity line centered at 121.83 MHz. The separation between the step, u c and the 
maximum, v a ^ increases at low temperature along with line broadening. 

Figure 9 shows the temperature dependence of K iso and K ax in NdCoPO in the 
temperature range 140 - 300 K (paramagnetic phase) along with those reported in 
LaCoPO [10] SmCoPO [T2] for comparison. With the decrease of temperature, both 
the shift and the line width increase due to the increase of magnetization. 

Inset of figure 9 shows the linear variation between Ki SO and K ax with \ which 
indicates unique isotropic and anisotropic hyperfine coupling constants. This further 
suggests that the presence of the impurity line in the NMR spectrum does not hamper 
the determination of the shift parameters of NdCoPO from theoretical fitting. Using 
equation 1 we obtain H™? and Hffi to be 7.14 kOe//i B and -2.34 kOe//i# respectively 
which are of same order of magnitude as in LaCoPO and SmCoPO shown in figure 
10. However, the 75 As NMR hyperfine coupling constant derived in case of NdFeAsO 
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Figure 9. (Color online) Temperature dependence of Ki SO and K ax . Filled circles, 
filled square, filled star corresponds to K iso for NdCoPO, LaCoPO [10] and SmCoPO 
[12] respectively. Open circles, open square, open star corresponds to K iso for 
NdCoPO, LaCoPO [10] and SmCoPO [12] respectively. Inset shows K iso and K ax 
versus M/H for NdCoPO. 



and NdFeAsOi-^F^ were one order of magnitude higher than that in NdCoPO [23]. 
This means the hyperfine coupling becomes weaker when As is replaced by P and Fe is 
replaced by Co. This can happen if such replacement alters the band structure so that 
the overlap of the different orbitals involved in H^f is affected. 

From the theoretical fitting of the K iso versus T curve using equation 3, the 
estimated value of P e jf obtained from the value of the Curie- Weiss constant (C) is 
2.23fi B with #=90 K for NdCoPO. For LaCoPO, P eff = IA/j b . This indicates that 
their is a contribution of Nd 4/ moment over that of Co-3d moment on shift even in 
the paramagnetic state. In case of NdFeAsOi_ x F x the P e ff value was close to local 
Nd-4/ moment where effect of Fe 3d is negligible [23]. However, in case of NdCoPO if 
we assume that the Co-3<i moment contribution is same as in LaCoPO, then from the 
Ration P«„ = ^ + (,/,,)» one can easily see that the P „ vaine to, Nd 4/ is 
lower than the local Nd 4/ moment, which reveals a signature of the itinerant character 
of Nd-4/ electrons in NdCoPO similar to that of the Co 3d electrons. 

Figure 10 shows that H™? follows linear relationship with P e ff from LaCoPO to 
SmCoPO. Whereas, H™? and Hffi do not decrease linearly with the decrease of lattice 
volume similar to that in As based series (where the 75 As NMR results are available only 
in LaCoAsO, CeCoAsO and the present results in PrCoAsO). In RECoPO series H % h s ? 
and Hffi decreases linearly with decrease of lattice volume from LaCoPO to NdCoPO. 
With further decrease of lattice volume in SmCoPO, both the coupling constants 
increase abruptly indicating a drastic change in the band structure. Therefore, the 
present results bring out the necessity of a systematic NMR study also in other members 
of i?i?CoAsO and RECoPO series, in order to get a clear microscopic picture about the 
magnetic properties of both the series. 
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Figure 10. (Color online) (a) Hffi versus P eff (fi B ) for LaCoPO [10], NdCoPO and 
SmCoPO [12 , (b) Hi 8 / (closed symbol) and H^ (open symbol) versus V (A 3 ). 



4-2. Nuclear spin-lattice relaxation rate (1/Ti) 
In general, {l/TiT) SF is given by 

{1/T X T) SF ex Y, \H hf (q)\ 2 X tt(q,tu n )/tu n 



(6) 



where x"((7> w n) is the imaginary part of the transverse dynamical electron spin 
susceptibility, 7 n and u n are the nuclear gyromagnetic ratio and Larmor frequency 
respectively. Hhf(q) is the hyperfme form factor. 

Nuclear spin lattice relaxation time (T"i) is measured at the high intensity peak 
(Pab in figure 8) of 31 P NMR spectrum at all the temperatures. Because of the presence 
of the impurity line near the position of the step (z/ c ) at almost all the temperatures, 
T\ is not measured at this position, to avoid any discrepancy in the measured values. 
The temperature dependence of (l/TiT) & in NdCoPO is shown in figure 11 along with 
those reported in LaCoPO [10], and SmCoPO [12] for comparison. It is interesting to 
note that the order of magnitude of relaxation time (T%) in LaCoPO and NdCoPO are 
same. Whereas, in SmCoPO it is two orders of magnitude shorter. This indicates that 
in case of NdCoPO, the contribution to T\ of the Nd-4/ electrons over that of Co 3d 
electrons, is not as strong as it is in SmCoPO, even if the effective moment of Nd-4/ 
is higher than that of Sm-4/. The abrupt enhancement of Hhf in SmCoPO (figure 
10) should be the reason for the dominant role of RE-4/ electrons in 31 P relaxation 
process of SmCoPO. So a comparison of the 31 P T\ data as well as the magnitude of 
the hyperfme coupling constants in LaCoPO, NdCoPO and SmCoPO clearly reveal a 
significant change in the electronic band structure in SmCoPO resulting an enhancement 
of H hf which shortens the spin-lattice relaxation time (Ti) along with the dominant 
effect of enhanced imaginary part of dynamic spin susceptibility. 

In the paramagnetic state, the P e ff value in NdCoPO is higher than that in 
SmCoPO still the value of T\ in NdCoPO is much longer that in SmCoPO. This should 
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be due to the decrement of the inter layer separation (csmCoPO < c^dCoPo), possibly 
affecting the band structure, (which could enhance the 31 P nuclear hyperfine coupling) 
and not due to the higher moment of Nd-4/ compared to that of Sm-4/ electrons. From 
the inset of figure 11, it can be easily seen that \jT{T at 300K in NdCoPO is higher than 
that in LaCoPO due to the enhanced P e // value, but for SmCoPO it becomes higher 
though P e ff value is small compared to NdCoPO. Moreover below the lattice volume of 
NdCoPO, the static magnetic properties (figure 10) as well as the dynamics of electrons 
changes drastically in case of SmCoPO. In the 1111 series Sm based systems always show 
interesting and different electromagnetic properties. In case of 1111 superconductors, 
SmFeAsOi-^Fz has highest T c and also the 19 F relaxation rate is order of magnitude 
higher than that of La or Pr based analogues due to the strong role of Sm 4/ electrons 
[241125]. The hybridization between Sm 4/ and Fe 3d is stronger than that in the other 
compounds and as a result the Sm 4/ electrons are more itinerant in nature. The same 
case also happens in Co based magnetic compounds. Further the order of magnitude of 
\jT\T in SmCoPO [12] and in SmFeAsOi_a;F x [251124] is same when probed at 31 P site 
and at 19 F site respectively, which indicates that the role of Co 3d or Fe 3d are almost 
suppressed by the dominant effects of Sm 4/ electrons. 

It would be interesting to probe the electromagnetic properties of the series 
Ndi-^Ya-CoPO, where increment of Y concentration may result in the decrease of the 
lattice parameters. It can happen that for a particular concentration of Y, the lattice 
parameters become same as SmCoPO. This study can expected to reveal directly that 
only lattice parameter (not the moment of rare earth) is responsible for the magnetic 
properties of these systems. 
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Figure 11. (Color online) Temperature dependences of (l/TiT) (, for LaCoPO (star) 
[IS], NdCoPO (circle) and SmCoPO (square) [12]. Inset shows \jT{T versus P eff for 
LaCoPO [TD], NdCoPO and SmCoPO [T2J. 



When the Knight shift and the nuclear spin-lattice relaxation process are governed 
by conduction electrons, \jT{TK 2 is constant. If there is an exchange interaction 
between the electrons then, using the Stoner approximation along with random phase 
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Figure 12. (Color online) T versus (l/T x TK) ab (half filled circle) and (\jT l TK' i l' 1 ) ah 
(filled circle) of NdCoPO. Inset shows (l/TiT) o6 versus K ab of NdCoPO. 



approximation, modified Korringa relation can be written as 
n(a), where So = [hj 47r/c s )(7 e / 'in) 2 and 

rv(a) = ((1 - a ) 2 /(l - « 9 ) 2 )f5- 



23(28] S /T 1 TK : 



2 
spin 



(7) 



a„ 



otQXoi. ) I 'x(O) is the g-dependent susceptibility enhancement, with a 



Xo(0)/x(0) representing the g = value. The symbol Qps means the average over 
q space on the Fermi surface. x(0) and Xo(?) represents the static susceptibility and the 
q mode of the generalized susceptibility of non-interacting electrons respectively. n(a) < 
1 means the spin-fluctuations are enhanced around q = 0, leading to the predominance 
of ferromagnetic correlations and n(a) > 1 signifies that spin-fluctuations are enhanced 
away from q = 0. This would indicate a tendency towards AF ordering (at q ^0). In 
case of LaCoPO k(o) < 1 both for ab plane and c direction, while for SmCoPO k{q) > 
1 for ab plane and K,(a) < 1 for c direction. In case of NdCoPO, n(a) a b < 1 (0.04). This 
indicates a ferromagnetic correlations in the ab plane. 




T/T„ 

c 



Figure 13. (Color online) P cS /P s vs T c /T plot (Rhodes- Wohlfarth plot). Solid line 
corresponds to Rhodes- Wohlfarth equation i.e. P c s/P s ex (Tc/To) -3 / 2 . The points 
(closed circle) are taken from [3B], LaCoPO from [TU], SmCoAsO from 37 j, SmCoPO 
from [12], LaCoAsO from [7] and NdCoPO: present work. 
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According to the SCR theory of weak itinerant ferromagnet, if 3D/2D spin- 
fluctuations are dominant [29[ 130] then \jT{T oc ^ 1 ( 3 / 2 ). Figure 12 shows the T versus 
{l/T\TK)ajb and (l/TiTK 3 ^ 2 )^ plots revealing dominant 3D FM spin-fluctuations in 
the ah plane of NdCoPO in the paramagnetic region in contrast to LaCoPO and 
SmCoPO where 2D FM correlations were present. It is to be pointed out that even if 
LaCoPO, NdCoPO and SmCoPO show a ferromagnetic transition for Co 3d electrons, 
the dimensionality of electron spin fluctuations of Co 3d are not same in all. In NdCoPO 
3D spin fluctuations of Co 3d dominates but it was 2D in case of LaCoPO and SmCoPO 
in the paramagnetic state. 

If 3D FM spin fluctuation dominates the relaxation process then [3H [H] 

(1/TiT) ~ 3h^ n H hf K/16^ B T . (8) 

where spin-fluctuation parameter T characterizes the width of the spin excitations 
spectrum in frequency space. Hence by plotting (l/TiT) a & versus K a b (inset of figure 
12), we have calculated T Q from its slope, which is 499 K, much smaller than that of 
LaCoPO and SmCoPO [12]. 

Further, according to the SCR theory of weak itinerant ferromagnet (WIF) 
[3 [321 [331 El] 

T c = (60c)- 3 / 4 P5 2 TJ /4 T 1/4 (9) 

where c=0.3353- • • and P s is M s0 in \i B unit. Ta which characterizes the energy width of 
the dynamical spin-fluctuation spectrum and evaluated to be 19730 K (using equation 
9) for NdCoPO. 

This SCR theory is also used to evaluate the coefficient F\ (an important spin 
fluctuation parameter) of M 4 term, in the Landau expansion of free energy, which can 
be written as [35] 



Fi = 4TJ/15T . (10) 

This equation indicates that there is a significant renormalization effects due to the 
zero point spin-fluctuations, which modify expansion coefficients of the free energy with 
respect to the uniform magnetization and further confirms that at T = OK, behavior 
of zero point spin-fluctuations have an important role in determining the magnetic 
properties in itinerant systems. F\ is 208028 K for NdCoPO (as estimated from equation 
10) which is greater than that of SmCoPO (57562 K), This suggests a more localized 
character of the d electrons in NdCoPO compared to that in SmCoPO. 

According to the SCR theory, in case of weakly itinerant ferromagnets, zero point 
spin-fluctuation mainly contributes to the total spin-fluctuations and the thermal spin- 
fluctuations increase with temperature above Tc, which makes the value of the ratio 
of effective paramagnetic moment to the saturation moment, P c g/P s greater than one. 
But in case of local ferromagnetic systems zero point spin-fluctuation is negligible and 
thermal spin-fluctuations are constant above Tc and which gives P e s/P s to be equal to 
one. Thus the ratio P e s/P s > 1 indicates that P c s does not contributes to the static 
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magnetic moment but to the dynamic magnetic response which indicates the role of spin- 
fluctuations in describing the magnetism of itinerant magnetic systems. From the above 
argument it is clear that the ratio P e s/P s is important to distinguish between itinerant 
and local systems. Further due to the presence of the spin-fluctuation in itinerant 
systems, the spin-fluctuation parameter T is also of interest and according to the SCR 
theory, Tq/Tq = 1 in case of localized ferro magnet, but as the itinerant character sets 
in, the value becomes less than one. It has been seen that all itinerant ferromagnetic 
materials follow a universal plot (Rhodes- Wohlfarth plot) between P e s/P s and Tc/Tq 
[35| [36] . To check the validity of the obtained spin fluctuation and thermodynamic 
parameters in NdCoPO and SmCoPO, we have plotted P e g/P s versus Tc/Tq in figure 
13. Though NdCoPO and SmCoPO are not purely itinerant ferromagnet and also 
there are contributions of 4/ electrons on Knight shift and relaxation rate (1/Xi) over 
Co 3d, we have tried to see whether NdCoPO and SmCoPO follow RW plot or not. 
From figure 13 one can see that NdCoPO and SmCoPO roughly follow RW relation 
P e fi/P s ~ (Tc/To) -3//2 . From the plot one can say that the derealization character of 
Co 3d electrons increases from LaCoPO to NdCoPO to SmCoPO. For comparison, we 
have also included SmCoAsO (37] in the RW plot, which is close to SmCoPO, indicating 
the change of the itinerant character of Co 3d is nominal for the replacement of As by 
P. 

5. CONCLUSION 

We have reported 75 As and 31 P NMR in PrCoAsO and NdCoPO respectively. The 
similarities and dissimilarities of the NMR results in the As and P based systems with 
different rare earths, are discussed. The role of 3d electrons, in the Co and Fe based 
systems, has also been revealed indicating an interesting interplay between 4/ and 3d 
in these magnetic systems. 

We have proposed a possible spin structure in PrCoAsO from 75 As NMR results. 
A scaling between vq and K iso or \ indicates the presence of a coupling between charge 
density and spin density fluctuations as proposed theoretically by Moriya for itinerant 
magnetic systems. It is shown that similar scaling also holds in the itinerant systems 
e.g. BaFe 2 As 2 [B], NaCo 2 4 [20] and LiFeAs [23]. 

The behavior of Knight shift and spin lattice relaxation data in the paramagnetic 
state of NdCoPO was compared with those reported in LaCoPO and SmCoPO. The 
role of RE-4/ effective moment is less prominent in NdCoPO. These results suggest 
a comparatively large modification of the band structure in SmCoPO, resulting from 
the decrement of inter layer separation. Nuclear spin lattice relaxation rate (1/Ti) in 
NdCoPO at 300 K, is nearly isotropic. The relaxation data was analyzed using the 
SCR theory of itinerant ferromagnet and was shown that, along the ab plane, three 
dimensional ferromagnetic correlations are present among the Co-3<i spins. Moreover, 
the different spin fluctuation parameters were evaluated. We conclude that the role of 
Co 3d electrons is not universal for different members of the family in the paramagnetic 
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state along the ab plane. 
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